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Abstract 

We consider two classical ensembles of the random matrix theory: the 
Wigner matrices and sample covariance matrices, and prove Central Limit 
Theorem for linear eigenvalue statistics under rather weak (comparing with 
results known before) conditions on the number of derivatives of the test 
functions and also on the number of the entries moments. Moreover, we 
develop a universal method which allows one to obtain automatically the 
bounds for the variance of differentiable test functions, if there is a bound 
for the variance of the trace of the resolvent of random matrix. The method 
is applicable not only to the Wigner and sample covariance matrices, but to 
any ensemble of random matrices. 

1 Introduction 

The Wigner Ensembles for real symmetric matrices is a family of n x n real 
symmetric matrices M of the form 

M = n-^/'^W, (1.1) 

where W = {^^^jj^^^-^ with w^^^ = w^^^ G K, 1 < j < < n, and 
^ 1^ j ^ k < n are independent random variables such that 

Here and below we denote E{.} the averaging with respect to all random param- 
eters of the problem. Let {Aj"^}"^^ - be eigenvalues of M. Since the pioneer 



work of Wigner [15] it is known that if we consider the linear eigenvalue statistic 
corresponding to any continuous test function 

n 

Afn[v]=Y.v{xf), (1.3) 
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then n ^J\fn[^p] converges in probability to the Hmit 

Hm n-^MnM = [ ip{X)Psc{X)dX. (1.4) 

n— ^-oo J 

where Psci^) is the famous semicircle density 

Psc{\) = ^\/4- A2l[_2,2]. 

The result of this type, which is the analog of the Law of Large Numbers of the 
classical probability theory, normally is the first step in studies of the eigenvalue 
distribution for any ensemble of random matrices. For the Wigner ensemble this 
result, obtained initially in [15] for Gaussian W = {wp,^}"^j^^-^^, was improved in 
|llj . where the convergence of Nn{X) to the semicircle law was shown under the 
minimal conditions on the distribution of W = {u^j^'' (^^^ Lindeberg type 
conditions). 

The second classical ensemble which we consider in the paper is a sample 
covariance matrix of the form 

M = n-^XX*, (1.5) 

where X is a n x m matrix whose entries 1^^^?^ ._, ,_, are independent 
random variables, satisfying the conditions 

E{^j?}=0, E{{X^l^f} = l. (1.6) 

Corresponding results on the convergence of normalized linear eigenvalue statis- 
tics to integrals with the Marchenko-Pastur distribution were obtained in [lOj. 

Central Limit Theorem (CLT) for fluctuations of linear eigenvalue statistics 
is a natural second step in studies of the eigenvalue distribution of any ensemble 
of random matrices. That is why there are a lot of papers, devoted to the proofs 
of CLT for different ensembles of random matrices (see [H El El El [U [T21 [HI E] ) . 
CLT for the traces of resolvents for the classical Wigner and sample covari- 
ance matrices was proved by Girko in 1975 (see [5] and references therein), but 
the expression for the variance found by him was rather complicated. A sim- 
ple expression for the covariance of the resolvent traces for the Wigner matrix 
in the case E{(wl^^)'^} = 2 was found in [8]. CLT for polynomial test func- 
tions for some generalizations of the Wigner and sample covariance matrices was 
proved in by using moment methods. CLT for real analytic test functions 
for the Wigner and sample covariance matrices was established in [2| under ad- 

,(")^2T. _ o 7T^r^„(n)N4\ _ q ^2 JC,,/") ^2^ 



ditional assumptions that ^{(< } = 2, E{{w)l')*} = 3E^{{w)l'y} = 3 (or 
= SEHiX^l^ 



^{(^!fc¥} = 3^M(^!r^)^} foi" the model In the recent paper [9] CLT 
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for the linear eigenvalue statistics of the Wigner and sample covariance matrix 
ensemble was proved under assumptions that E{{wl^^ )'^} = 2, the third and the 
forth moments of all entries are the same, but E{{w^j^)^} is not necessary 3. 
Moreover, the test functions, studied in [[DJ, are not supposed to be real analytic. 
It was assumed that the Fourier transform ip of the test function ip satisfies the 
inequality 

j{l + \kf)\!p{k)\dk<oo, (1.7) 

which means that (p has more than 5 bounded derivatives. 

In the present paper we prove CLT for the Wigner ensemble ()l.ip under the 
following assumptions on the matrix entries 

^{i^tkf] = "^4, sup sup E{\w^;i^\^+''} = Wi+e, <oo, ei> 0. (1.8) 
n l<j<k<n 

We consider the test functions from the space Tig, possessing the norm (cf ()1.7p ) 
11^11"^ = Jil + llkinmi'dk, s>3/2, m = ^Je"'^v{x)dx. (1.9) 



Theorem 1 Consider the Wigner model with entries satisfying condition U.8\) . 
Let the real valued test function ip satisfy condition ||v'||3/2+e < oo (e > 0). Then 
A/'°[(/3] converges in distribution to the Gaussian random variable with zero mean 
and the variance 




where K4 = — 3. 

Let us note that similarly to the result of [9j it is easy to check that Theorem [T] 
remains valid if the second condition of (jl.Sp is replaced by the Lindeberg type 
condition for the fourth moments of entries of W 

lim l'^Ht) = 0, Vr > 0, (1.11) 

where 

1 " 

j,k=i ' 
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The proof will be the same as for Theorem [H but everywhere below n"^^/^ will 
be replaced by L„(t)/t''', with some positive 7. 

The proof of Theorem [1] is based on some combination of the resolvent ap- 
proach with martingal bounds for the variance of the resolvent traces, used before 
by many authors, in particularly, by Girko (see [S] and references therein). An 
important advantage of our approach is that it is shown by the marginal difference 
method that (see Proposition [2] below) 

Var{TrG(z)} < C7/|9z|^ G{z) = {M - z)-\ (1.13) 

while in the previous papers the martingal method was used only to obtain the 
bounds of the type Var{Tr G{z)^^} < nC{z). The bound (|1.13p will be combined 
with the inequality 

Proposition 1 For any s > and any M 

Var{Mn[ip]} < CsWifWi dye-yy^'-^ \ar{TY G{x + iy)]dx. (1.14) 

Jo J~OD 

The proposition allows one to transform the bounds for the variances of the 
resolvent traces into the bounds for the variances of linear eigenvalue statistics 
oi (p € Tis, where the value of s depends on the exponent of \Qz\ in the r.h.s. of 
(jl.lSp . It is important, that Proposition [T] has a rather general form and therefore 
it is applicable to any ensemble of random matrices for which the bounds of the 
type ()1.13p (may be with a different exponent of I'^zl) are found. This makes 
Proposition [1] an important tool of the proof of CLT for linear eigenvalue statistics 
for different random matrices. The idea of Proposition [T] was taken from the paper 
[7], where a similar argument was used to study the first order correction terms 
of n~^E{A/'n [</']} for the matrix models. Having in mind Proposition [H one can 
prove CLT for any dense in Tis set of the test functions, and then extend this 
result to the whole Tis by the standard procedure (see Proposition [3]) . In the 
present paper for this aim we use a set of convolutions of integrable functions 
with the Poisson kernel (see (I2.32P and (12. 3p ). This choice simplifies considerably 
the argument in the proof of CLT and makes the proof more short than that in 
the previous papers [Il[2l[9]. 

The result for sample covariance matrices is very similar. We assume that 
the moments of the entries of X from (jl.Sp satisfy the bounds 

E{(X|;;))4} = X4, sup sup E{|xg)|4+^i} = X4+ei < 00, ei > 0. (1.15) 

n l<j<k<n 

Theorem 2 Consider a random matrix \1.5\) - /il.6\) with entries of X, satis- 
fying the condition \1.15\) . Let the real valued test function satisfy condition 
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llvll3/2+e < oo (e > 0). Then M^i'p] in the limit m,n — )• oo, m/n — )■ c > 1 
converges in distribution to the Gaussian random variable with zero mean and 
the variance 

1 ra+ f"-+ / Aip\'^ (Ac- {Xi - am){X2 - am)]dXidX2 



2vr2 7,_ 7„_ VAA; ^Ac - (Ai - a^) V4c - (A2 - a^y 

where — = — t-^—^, 1^4 = X4 — 3 is the fourth cumulant of entries of X, 
AA Ai — A2 

a± = (1 ± \/c)^, and Om = + a ). 

2 Proofs 

Proof of Proposition [H Consider the operator Vs 

vrf{k) = {i+2\k\ym (2.1) 

It is easy to see that for fixed n Var{A/'n[<^]} is a bounded quadratic form in the 
Hilbert space T-L of the functions with the inner product {u,v)s = (VsUjVsv), 
where the symbol (., .) means the standard inner product of L2(M). Hence there 
exists a positive self adjoint operator V such that 

Var{AA„[(^]} = (Vif,^) = TV(n^Vn^) 

where 11^ is the projection on the vector ip 

{U^f){x) = ip{x){f,^)\M^' 

where ||.||o means the norm ()1.9p with s = 0. We can write 



But it is easy to see that 

(p,n^/)(x) = (i?,^)(x)(/,(^)||(^|loS 

hence, 

||Psn<^|| = ll^^syjllo = WvWs- 

Therefore we can write 

Var{AA„M} = TriU^VsV^^VV^^V.U^) < \ \VsU^\\^Tr(V;^VV;^) (2.2) 
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But since for any u,v £ L2(M) we have 

r{2s){V-^u,v) = r(2s) J {1 + 2\k\)-'^'u{k)^k}dk 



oo 



dye yy 







-^\^\yu{k)v{k)dk = / dye-yy^'-\Py*u,Py*v) 
Jo 



OO err 

dye-yy^'~^ dx Pyix - \)Py{x - fi)u{X)lJ(jL)d\dn, 



where the symbol * means the convolution of functions, and Py is the Poisson 
kernel 

^ -K^x^ + y ) 

This implies 

T{2s)V-\\, fi) = dye-yy^'~^ j dxPy{x - X)Py{x - /z), (2.4) 

and so 

T{2s)1:t{V-^VV-^) = dye~yy^'~^ j dx(vPy{x - Py{x - .)) 

= dye-yy^'-^ J dxVar{Mn[Py{x - .)]} 

= J dye-yy'^'-^ j da;Var{$>Tr G(x + %)}. 

This relation combined with (|2.2p proves (|1.14p .D 

In what wallows we need to estimate E||t(;j^^|^} (see the proof of Proposition 
[2]). Hence, if £i < 4, then it is convenient to consider the truncated matrix 

mW = {m;W}^^.^^, = m,,1|a,,^,|<„ mM°-e{mW}. (2.5) 



Lemma 1 Let Mri[^\ = Tr (^(M^"^)") he the linear eigenvalue statistic of the ma- 
trix M^'^^° , corresponding to the test function (p with bounded first derivative. 
Then 
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Proof. Consider the matrix M{t) = M + t(M — M). Let {\i{t)} be eigenvalues 
of M{t) and be corresponding eigenvectors. Then 

pi 

< WWooJ^ dtE{^|(M'(t)V^,(t),Vi(t))|} < llv^'lUEjTr |M-M|} 
= 1 1 1 1 oo e{ I ^ u^^. (M - M)ijujk I } 

k ij 

< ||(/p'||ooE{ ^ |(M - M)ij-|} <sup|y.'|L„ (r)/r3, 

ij 

where M'(t) = ^M(t) = (M - M), U = {uik} is the unitary matrix such that 
M - M = U*AU, where A is a diagonal matrix and \M - M\ = U*\A\U. Hence, 

|gixAr°M _ gixAr„M°| < 2Pr{M(^) ^M} + |x| (e{AA„ [(/?]} - E{7V;[v9]}) 
<o(1) + C7|x|||v;'||ooL„(t)/t^ 

□ 

It follows from Lemma [T] that for our purposes it suffices to prove CLT for 
Af°[ip]. Hence, starting from this point we will assume that M is replaced by 
M^"^^", but to simplify notations we will write M instead of M^'^^° just assuming 
below that the matrix entries of W satisfy conditions 

B{w,k} = 0, B{w%} = 1 + o(l), (j / k), B{w]^} =W2 + o(l), (2.6) 
B{w%} = W4 + 0{1), 

B{\wjkf} < W4+ey-'''\ B{\Wjk\^] < W^+ey-''/^. (2.7) 

Here and below we omit also the super index (n) of matrix entries Wj'^ and Xj^^ . 

Proposition 2 // the conditions \2. 6]) are satisfied, then for any 1 > 5 > 
n 

Var{7„} < Cn-i J^E{|Q,(z)|^+'^}/|9z|3+^ Var{7„} < C/\'^z\\ (2.8) 

i=l 

If the conditions of \2. are also satisfied, then 

E{|7°|^} < Cn-i-^i/Vl^-^l^^- (2-9) 
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Proof. Denote E<fc the averaging with respect to {t(^ij}i<i<j<A:- Then, according 
to the standard martingal method (see [3]), we have 



Var{74 = ^E{|E<fc_i{7„} -E<fc{7„}|2}. 



(2.10) 



k=l 



Denote E^, the averaging with respect to {wki}i<i<n- Then, using the Schwarz 
inequality, we obtain that 

|E<fc_i{7n} - E<fc{7„}|2 = |E<fc_i{7n - Efc{7n}| < E<fc_i{|7„ - Efc{7„}|2}. 

Hence 



Var{7„}<^E{|7„-Efc{7„}|2}. 



(2.11) 



k=l 



Let us estimate the first summand (with k = 1) of the above sum. The other 
ones can be estimated similarly. Denote M^^^ the (n — 1) x (n — 1) matrix which 
is the main bottom (n — 1) x (n — 1) minor of M 



Gd) = (M(i) - z)-i 
We will use the identities 

TrG-TrG«^ l + (G«G«m«,m« 



,«;in)GM"-^ 

1 + B{z) 



G 
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-A 



-1 



z + n-^/'^wii + (G(i)m(i), m(i)) 



Aiz) 



(2.12) 



(2.13) 



where (., .) means the standard inner product in ^. 

The first identity of (j2.13p yields that it suffices to estimate E{|i?j4^^ — 
Ei{SA"^}p} and E{|A"^ - Ei{^"^}p}. We wih estimate the first expression. 
The second one can be estimated similarly. Denote = ^— Ei{^} for any random 
variable ^ and note that for any a independent of {wu we have 

E^m'}<E^{\C-a\'}. 



Hence it suffices to estimate 

B Ei{B} _ B° 
A ~ B,{A} 



Al B 



< 



Bl 



Ei{A} 



+ 



Al 



9zEi{^} 



Ei{A} Ei{^} A 
Let us use also the identities that follow from the spectral theorem 

9(GWmW,m(i)) = 9z(|G(i)| V^^m^^)), QTr G^^^ = QzTr jG^^^I^. (2.14) 
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where = (G^^^G^^)*)^/^. The first relation yields, in particular, that |B/yl| < 

Moreover, using the second identity of (j2.14p . we have 



n 



)|2 



|z + n-iTr G(i)|2 
Since 



(n-^Tr |GWp)^(n-iTr |GWp)i-^ 
+ n-iTr G{i)|i+'5|z + n-iTr GWp- 



< G- 



Ei{A}|i+5- 
(2.15) 



(2.16) 



Ei{|^?|2} < Gn-^Tr jG^p + Gn'S 
we get by (j2.15p and the second identity of ()2.14p : 



-l-<5 



Similarly 



El 



Ei{A} 



Gn-2Tr|G«|4 Gl^-i-^--^ 

< ^ < - 



z| -n -Tr |G(i)| 



+ n-iTr G(i)|2 



+ n-iTr G{i)|2 



< G 



(2.17) 



|Ei{A}| 



l+<5 



Then, using the Jensen inequality |Ei{^}| ^ < Ei{|A| ^}, and the second iden- 
tity of (j2.13p . we conclude that 



E{|(7„(.))?|^}< 



G 



E{|Gn(.)|^+n- 



Then (j2.1ip implies 

To prove ()2.9p we use the inequality similar to (|2.1ip (see [1]) 

n 

E{|7°in<C^^EE{l7n-Efe{7n}r}- 

k=l 

Thus, in view of (|2.13p . it is enough to check that 

Ei{|^?|^} < Gn-l~"l/2|c^z^^ ^i{\B°\\} < Gn-i-^i/2|cj_^| 



(2.18) 



(2.19) 



The first relation here evidently follow from (j2.16p , if we take the forth degree of 
the r.h.s., average with respect to {wu}, and take into account (j2.7p . The second 
relation can be obtained similarly. □ 

Proposition [2] gives the bound for the variance of the linear eigenvalue statis- 
tics for the functions 93(A) = (A — z)~^ . We are going to extend the bound for a 
wider class of test functions. 
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Lemma 2 // ||</'||3/2+e ^ oo, with any e > 0, then 

Vsir{Mn[ip]}<C,\Ml/2+s- (2-20) 
Proof. In view of Proposition [1] we need to estimate 

/oo 
Var{7„,(x + iy)}dx 
-oo 

Take in (|2.8p 5 = e/2. Then we need to estimate 

/oo 
E{|G,j(x + iy)|i+^/2|^x, j = l,...,n. 
-oo 

We do this for j = 1. For other j the estimates are the same. The spectral 
representation 

Nn{dX) 



Gn 

and the Jensen inequality yield 



X — X — iy 



Taking s = 3/2 + e in (fTTiD . we get 

/>oo 

□ 

To simplify formulas we will assume below that {wjk}i<j<:k<n are i.i.d. and 
{'Wjj}i<j<n are i.i.d. Note that this assumption does not change the proof seri- 
ously, it just allows us to write the bounds only for Gn instead of all Gu. 

The next lemma collects relations which we need to prove CLT. 

Lemma 3 Using notations of (2.13\) we have uniformly in zi,Z2 ■ 9-Zi,2 > cl with 
any a > 0: 

E{{A°f}, E{|^°|^}, E{{B°f}, E{\B°\^} = 0{n-^-''/^), (2.21) 
Bi{A°{zi)A°{z2)} = -TrGW(zi)GW(z2) + W2 (2.22) 



n 



n 



i 

nEi{A^zi)B°{z2)} = n-^Ei{A°{zi)A°{z2)}, (2.23) 

dZ2 
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Var{nEi{yl°(zi)^°(z2)}} = O(n-i), (2.24) 
Var{riEi{A°(zi)S°(z2)}} = 0{n-^), 

n\ln\z) - Inizt} = 0{n-'-'^/^). (2.25) 

Moreover, 

Var{Gi^)(zi)} = 0{n-'), \E{G^ (z^)} - E{G..(zi)}| = 0{n-') (2.26) 

|E{7i')(z)}/n - f{z)\ = 0{n^'), \B-'{A{z)} + f{z)\ = 0{n''). (2.27) 

Proof. Note that since ^^^(G^^^m, m) > 0, we can use the bound 

> ^ \A'^\ < < (2.28) 
Relations (|2.21|) fohow from the representations 

^o^^o^^-io(i)(^)^ BO = Bl + n-^^^t\^), (2.29) 

combined with ([2T9]) . and ([2:9]) . apphed to Relations ([2:22]) and ([2:23]) 

follow from (|2.16p and (|2.29p . if we take the products of the r.h.s. of (|2.16p with 
different z and average with respect to {wu}. Relation ([2.25p follows from ([2.13p . 
([Ml]), and ([2:28]) . The first relation of is the analog of the relation 

Var{Gii(zi)} = Var{Gii(zi)} = ©(n^^) (2.30) 

if in the latter we replace the matrix M by M^^\ But since Gii(zi) = —A~^{zi), 
dOO]) follows from ^ni\\ and (ICT]1 . The second relation of (ITM]) follows from 

The first relations of ([2.27p follows from the above bound for n^^E{'~fn — In^} 
and the well known estimate (see e.g. [8j) 

n-iE{7„} - fiz) = O(n-i). 

The second one of (|2.27p is the corollary of the above estimate and of the relation 

E-H^(^)} = {z + E{7«}/n)-i = {z + f{z))-' + 0{n-') = -f{z) + 0{n-'). 

Finally we obtain the first bound of ([2:2^]) from ([2:22]) . ([2:261) . ([125]), and the 
identity 

T..OW(..)Oa)fe) = TV°"'<")-°"'''-\ (2^31) 

Zl - Z2 

The second bound of p.24|) follows from the first one, ([2.23p . and the Cauchy 
theorem. □ 
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Proof of Theorem [IJ We prove first Theorem [T] for the functions iprj of the 
form 

f-n = Pri*^o, J \vo{^)\d\ < C < oo, (2.32) 
where Prj is the Poisson kernel defined in ()2.3p . One can see easily that 

^nif] = - I Mf^)^lnizt^)dfJ', z^ = fi + irj. (2.33) 



Set 

Zn{x) = Eje^^^-t"^]}, e{x) = e^^^-t"^!, Yn{z,x) = B{TtG{z)e° (x)}. (2.34) 
Then 

d X f 

— Z„(x) = ^ / mi^J-)(Xiz^^,x) - Y{z^,x))dn. (2.35) 

On the other hand, using the symmetry of the problem and the notations of 
()2.13p . we have 



Yn{z, x) = E{Tr G{z)e°{x)} = n'E{Gu{z)e° (x)} (2.36) 
= -nE{(^-i)°ei(x)} - nE{{A~'ne{x) - ei(x))} = T, + T2, 

where 

eiix) = e-(^"-i[^l)°, (Mi'lMT = (Tr9.(M(i)))° = J d^M^)'iln\z,). 
Let us use the representation 

E{A} E2{yl}^E3{A} E4{^} ^ ^E4{^} " ^ ' 

Since ei(x) does not depend on {wu}, using that E{...} = E{Ei{...}}, we obtain 
in view of (12371) and (fi:2Tl) 

_ E{nEi{^°(z))}e^(x)} E{nEi{(A°(z))2}e^ (x)} , 
^" E2{^} E3{A} 

Relations (|2.24p implies 

|E{nEi{(^°(z))2}e°(x)}| < Vari/2{nEi{(^°(z))2}}Vari/2{e°(x)} = 0{n-^/^), 
thus 

Ti = E-2{^}E{7«e^(x)} + 0(n-^i/2). 



12 



But the Schwarz inequality and (|2.25p yield 

|E{(7i'))°ei(x)} - E{7°e(x)}| < E{|(7«)° - 7;^|(1 + |x||7°|)} (2.38) 
< C{1 + |x|)Ei/2{|(7«)o _ ^o|2| ^ o{n-'l^). 

Thus, we have 

Ti = E-2{A(z)}y„(z, x) + 0{n-^'^) = f{z)Yn{z, x) + 0{n-'-''^). (2.39) 
To find T2, we write 

e(x)-ei(x) = zx j ^o(/.)(9(7° - 7n^)ei(x) + ©((7^ - (7!'^)°)'))^/^. 

Using the Schwartz inequality, ()2.13p . (|2.25p . and (|2.22p . we conclude that the 
term 0((7^ — (7!^^)°)^) gives the contribution 0{n~^^/^). Then, since ei(x) does 
not depend on {wu}, we average first with respect to {wu} and obtain in view 
of ([2:2^ 



^2 = -^y (i/i¥'o(^)E{ei(x)(^-i)°(z)9(7„-7i'T(^^)} + 0(n-^^/') 
= d^^o(/.)E{ei(x)Ei{(^-i)°(.)9(i±^)°}}+0(n--/^) 

= dM<^o(/.)E{(^-i)°(.)s(i±^)°}E{ei(x)} + 0(n--/^). 

Using (I2.37P and (I2.2ip . we conclude that only linear terms with respect to B° 
and A° give non vanishing contribution, hence in view of (j2.23p and ()2.24p we 
obtain 

Dn{z,z^) ■.= n^[{A~^r{z){{l + B{z^))A-\z^)y] 

^ (l + E{i?(z^)})nEi{A°(z)^°(z^)} _ n^i{A°{z)B°{z^)} 

W{A{zm^{A{z,)} ^^{A{z)}^{A{z,)}^ ^ 

= /2(z)/2(z^)(l + /'(z^))(E|-TrG«(zi 



i 

+ ?E{5:G,«(z)G«(.,)})+0(n--/^), 
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where we used also (^27\) to replace E-^{A{z)} by f{z) and E{B{Zf,)} by f'{Zf,). 
The identity (fOT]) yields 

Dr^{z,z^) :=2/2(z)/2(^^)(i + /'(^^))(M^/M + ^2/2) (2.40) 

V Z Z^ / 

+ 2/^(=)/(.,)A(M^.)- 

dz,, V 2; - z,, 



+ ^^A{f{z)f\z,){l + /'(z^)) + f\z)f{z,)f{z^)) + 0(n--^/2). 

In addition, similarly to (I2.38p . we have 

E{ei(x)} = Z,(x) + 0(n-i/2). 
Hence, relations (f236D - (inoD imply 



Yn{z,x)=f {z)Yn{z,x) +ixZn{x) J dflip{n) + o(l), 

y„(z,x) =ixZ„(x) I (i/.v.oM ^"^^'^^^^T^^"^^'^^ +o(l), (2.41) 



Using the relations 



/(.)(/'(.) + l) = :r^ = -^. /'- 



1-/2(2) v/F^' ' v/F^' 

we can transform C„(z, z^) to the form 



''^ (Z-Z^)2\^(^2_ 4)1/2(^2 _ 4)1/2 y ^ (,^2 -4)1/2(^2 _ 4)1/2 

Now, taking into account (I2.35p . (I2.4ip . and ()2.42p . we obtain the equation 

-^Znix) = -xV[^o, v]Zn{x) + o(l) (2.43) 
ax 

^[¥^0,??] = ^ j j (po(w)v'o(^2)(c(z^,,z;:7) + C(z;i7,z^2) 

- c{z^^,z^^) - c(z;rr,^))f^/"if^^2- 



14 



Now if we consider 

then ()2.43p yields for any |x| < C 

= o(l), 

and since Zn{0) = Z„(0) = 1, we obtain uniformly in |x| < C 

Zn{x) = 1 + 0(1) 

= e-^-'^l^'']/^ + o(l). (2.44) 

Thus, we have proved CLT for the functions of the form ()2.32p . To extend CLT 
to a wider class of functions we use 

Proposition 3 Let {^j"^}"^]^ be a triangular array of random variables, Mn[^] = 

n 

^^93(^1'"'^) be its linear statistics, corresponding to a test function : M — t- M, 

1=1 
and 

Vn[^] = Var{AA„[^]} 

be the variance ofJ\fn[^p]- Assume that 

(a) there exists a vector space C endowed with a norm ||...|| and such that Vn 
is defined on JC and admits the bound 

Vniif] < C\\ip\\\ yep e C, (2.45) 

where C does not depend on n; 

(b) there exists a dense linear manifold Ci <Z C such that CLT is valid for 

J\fn[f], f £ Ci, i.e., if Zn[xip] = E|e*^-'^"['^] | is the characteristic function of 

o 

n~^''^J\fn[^], then there exists a continuous quadratic functional V : L\ ^ M+ 
such that we have uniformly in x, varying on any compact interval 

lim = V99 G £i; (2.46) 

n— ^-oo 

Then V admits a continuous extension to C and CLT is valid for all Nn['-p\, 
Proof. Let {(fk} be a sequence of elements of £1 converging to ip £ C We 

o 

have then in view of the inequality le^"" — e^^\ < \a — b\, the linearity oi J\fn[^] in 



15 



the Schwarz inequality, and ()2.45p : 



Znixf) - Zn{xip)\^=^^ 



1^ u 



< |x|Var^/^{7V„[(/7 - ipk]} < C\x\ \\ip - ipk\\. 



Now, passing first to the hmit n — )• oo and then k — t- oo, we obtain the assertion. 

□ 

The proposition and Lemma [2] allow us to complete the proof of Theorem [TJn 

Proof of Theorem The proof of Theorem [2] can be performed by the same 
way as that for Theorem [TJ We start from the proposition which is the analog of 
Proposition [21 

Proposition 4 Let 7„ = Tr G{z), where G{z) = {M — z)~^ and M is a sample 
covariance matrix il.5\} with entries satisfying U.6\) and U.15\) . Then inequalities 
(TU) hold. 



Taking into account Proposition U on the basis of Proposition [T] and Lemma 
[2] we obtain immediately the bound (I2.20p for the variance of linear eigenvalue 
statistics of sample covariance matrices. Then one can use the same method as 
in the proof of Theorem [T] to prove CLT for 99^ of (I2.32p or just use the result 
of [9] for the functions, satisfying conditions (jl.7p . Then Proposition [3] implies 
immediately the assertion of Theorem [2j 

Thus, to complete the proof of Theorem [2] we are left to prove Proposition [H 

Proof of Proposition^ Similarly to the proof of Proposition [2] we use the iden- 
tity (j2.10p where this time E<jt means the averaging with respect to {^j7}i=i,.,m,j<fc- 
Then we obtain (j2.1ip with meaning the averaging with respect to {-^fcz};=i,.,m- 

Denote M^^-* = X^^'^ X^^'^* , where the (n — 1) x m matrix X^^^ is made from 
the lines X, from the second to the last one. Then denote 

= (m(i) - z)-\ = (M12, . . . Min) 

and use (I2.13P with these G^^^ and m^^\ 

To obtain the estimate for Ei||7„ — Ei{7„}p| we need (as in the proof 
of Proposition E]) to estimate Ei{| p}/(9zEi{^})2 and Ei{\Bl\^} / (Ei{A}f . 
Since G^^^ does not depend on {Xij}j=i^.^m, averaging with respect to {^ii}i=i,.,m) 
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using the Jensen inequality and (|2.47p . we get 

Ei{^} = ilYG^MW, 9Ei{^} = — TYG(i)mWg«*, 
n n 

< ^Ei-^{n-iTrG«M(i)G«*}Ef{n-iTYG(i)(M(i))(i+^)/^G«*} 

< _^Ej-'5|n-iTYG(i)M«G«nEf|n-iTr(M«)(i+^)/n. 

But it is known (see [3j and references therein) that for any fixed 5 > 

E{n-iTr(M«)(^+'^)/^} < {2 + + 0(1). (2.47) 

Combining this bound with the above inequahty and repeating the argument of 
Proposition H we obtain the bound j^HT^i . The bound for Bi{\B^\'^} /BHA} can 
be obtained similarly. 
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